Combined GW and dynamical mean field theory: 
Dynamical screening effects in transition metal oxides 
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We present the first dynamical implementation of the combined GW and dynamical mean field 
scheme ("GW+DMFT") for first principles calculations of the electronic properties of correlated 
materials. The application to the ternary transition metal oxide SrVOs demonstrates that this 
schemes inherits the virtues of its two parent theories: a good description of the local low energy 
correlation physics encoded in a renormalized quasi-particle band structure, spectral weight transfer 
to Hubbard bands, and the physics of screening driven by long-range Coulomb interactions. Our 
data is in good agreement with available photoemission and inverse photoemission spectra; our 
analysis leads to a reinterpretation of the commonly accepted "three-peak structure" as originating 
from orbital effects rather than from the electron addition peak within the t2g manifold. 
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Describing electronic excitations of materials with 
strong electronic Coulomb interactions remains one of the 
challenges of modern condensed matter theory [l| . Pho- 
toemission satehite features, quasi-particle mass renor- 
malizations or - in the strong coupling limit - the local- 
ization of electrons are phenomena well beyond band the- 
ory. Within effective low-energy models, such as the Hub- 
bard or Anderson lattice model, insight has been gained 
into the low-energy physics of correlated materials on a 
qualitative level. However, when relating this knowledge 
to specific materials, attempting a quantitative descrip- 
tion, several challenges persist: The first one is a question 
of energy scales: in fact, much of the physics of interest 
in correlated materials takes place on the scale of several 
eV, well beyond what is captured within a low-energy de- 
scription. Optical properties of transition metal oxides, 
to cite a specific example, are determined by the onset 
of optical transitions involving ligand states already at 
scales of only a few eV. A second issue is the determi- 
nation of the model parameters, e.g., the effective low- 
energy hopping parameters and Coulomb interactions, in 
particular the local "Hubbard t7". These two issues are 
in fact intimately related, since high energy screening 
processes renormahze these parameters, and an appro- 
priately downfolded low-energy model should take this 
into account 0-01 • 

In this Letter, we present the first dynamical imple- 
mentation of a combined GW and dynamical mean field 
( "GW-I-DMFT" ) scheme that addresses the above issues, 
by treating explicitly both, dynamical Coulomb interac- 
tions on the "correlated" manifold of orbitals, and correc- 
tions to the band structure of degrees of freedom that five 
on higher energy scales. The application to a well-studied 



correlated metal, SrVOa, shows some expected features: 
a renormahzation of the quasiparticle band structure, the 
formation of a lower Hubbard band and the correction of 
the position of the oxygen p-states. It also leads to less 
expected features, namely the absence of a clearly distin- 
guishable upper Hubbard band in the total spectral func- 
tion along with a reinterpretation of the feature seen at 
about 2.7 eV in inverse photoemission that we associate 
with the Qg states. This attribution is in agreement with 
cluster model calculations Q and allows to resolve appar- 
ent contradictions between cluster and dynamical mean 
field methods for this compound. In the orbital-resolved 
spectral function we furthermore predict additional satel- 
lite features induced by the frequency-dependence of the 
effective Coulomb interaction. 

The combined GW-I-DMFT scheme was proposed a 
few years ago, in order to avoid the ad hoc nature of 
the Hubbard parameter and the double counting inher- 
ent to conventional combinations of dynamical mean field 
theory with density functional theory in the local den- 
sity approximation (so-called "LDA-I-DMFT" schemes 0, 
llOll). The starting point is Hedin's GW approxima- 
tion (GWA)[lH. in which the self-energy of a quan- 
tum many-body system is obtained as a product of the 
Green's function G and the screened Coulomb interac- 
tion W = e^^V. The dielectric function e, which screens 
the bare Coulomb potential V, is -within a pure GW 
scheme- obtained from the random phase approxima- 
tion. The GW-I-DMFT scheme, as proposed in jlj], com- 
bines the first principles description of screening inher- 
ent in GW methods with the non-perturbative nature of 
DMFT, where local quantities such as the local Green's 
function are calculated to all orders in the interaction 
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from an effective reference system ("impurity model"). 
In DMFT, one imposes a self-consistency condition for 
the one-particle Green's function, namely, that its on- 
site projection equals the impurity Green's function. In 
GW-I-DMFT, the self-consistency requirement is general- 
ized to encompass two-particle quantities as well, namely, 
the local projection of the screened interaction is required 
to equal the impurity screened interaction. This in princi- 
ple promotes the Hubbard U from an adjustable param- 
eter in DMFT techniques to a self-consistent auxiliary 
function that incorporates long-range screening effects 
in an ab initio fashion. Moreover, the theory provides 
momentum dependence to quantities (such as the self- 
energy) that are local within pure DMFT. Despite these 
conceptually attractive features, however, no dynamical 
implementation, even at the non-selfconsistent level has 
been achieved so far (see however the static implemen- 



12l |. as well as model studies of a one-band 



tation in ^ 

extended Hubbard model in 13, 1^). 

From a formal point of view, the GW-f DMFT method, 
as introduced in Ref. [12[, corresponds to a specific ap- 
proximation to the correlation part of the free energy of 
a solid, expressed as a functional of the Green's function 
G and the screened Coulomb interaction W : the non- 
local part is taken to be the first order term in W, while 
the local part is calculated from an impurity model as in 
(extended) dynamical mean field theory. This leads to a 
set of self-consistent equations for G, W, the self-energy 
S and the polarization P [l5|, |l6[ . Specifically, the self- 

where the 
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energy is obtained as S = '^locai + 
local part Tiiocai is derived from the impurity model[56j. 
In practice, however, s- or p-orbitals are rather delocal- 
ized. Here, we therefore propose a practical scheme, in 
which only the local part of the self-energy of the "cor- 
related" orbitals is obtained from an impurity reference, 
and all other local and non-local components are approx- 
imated by their first order expressions in W. Green's 
functions and self-energies are expressed as matrices in a 
basis of atom-centered Wannier functions 57[ . For prac- 
tical reasons, we avoid the costly GW self-consistency 
cycle, and in particular the update of the dynamical in- 
teraction U{uj): Instead, we calculate U{bj) from the con- 
strained random phase approximation (cRPA) 0, 
and keep it fixed when iterating the local one-particle 
quantities to self-consistency. [58 1 

Our target material, SrVOa, crystallizes in the cubic 
perovskite structure: the V^^ ions are surrounded by 
oxygen octahedra, occupying the sites of a simple cubic 
lattice. The Sr^+ cation sits in the center of the cubes. 
The electron count leaves a single d electron in the V- 
d states, which is largely responsible for the electronic 
properties of the compound. The octahedral crystal field 
splits the V-d states into a lower-lying threefold degener- 
ate t2g manifold, thus filled with one electron per V, and 
an empty Cg doublet. The compound exhibits metallic re- 
sistivity with a Fermi liquid T^ behavior up to room tem- 



perature [l9| and temperature-independent Pauli param- 
agnetism without any sign of magnetic ordering |20i] . Hall 
data and NMR measurements [18,l21| as well as angle- 
resolved photoemission spectra [22| that find the Fermi 
surface of SrVOa as cylindrical Fermi sheets, in agree- 
ment with theory, confirm the picture of a normal Fermi 
liquid. 

These properties make SrVOa an ideal model material 
for studying the effects of electronic Coulomb interactions 
[il [23I- I35I . Recent works include, e.g., studies of surface 
effects |36j . kink structures in the momentum-resolved 
spectral function 37-39|, or non-local correlation effects 



401 |4l| in relation to laser ARPES [42] which finds a 



"dip" at the Fermi level or a maximum of the quasi- 
particle peak slightly below (at around -0.2 eV). The ex- 
perimentally observed quasiparticle weight is about 0.5- 
0.6 [isl. I43l - l49| ■ with the remaining spectral weight being 
carried by a photoemission (PES) (Hubbard-) satellite at 
around -1.6 eV 47, 50]. Inverse photoemission (BIS) has 
located the electron addition d^ peak at an energy 

of about 2.7 eV 

Figure [1^ summarizes the LDA electronic structure: 
the 0-p states disperse between -2 and -7 eV, separated 
from the t2g states whose bandwidth extends from -1 eV 
to 1.5 eV. While the t2g and Cg bands are well separated 
at every given k-point, the partial density of states (DOS) 
slightly overlap, and the eg states display a pronounced 
peak at 2.3 eV. Finally, peaks stemming from the Sr-d 
states are located at 6.1 eV and 7.1 eV. We have super- 
imposed to the LDA DOS the experimental PES and BIS 
curves taken from 
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The comparison reveals the 
main effects of electronic correlations in this material: as 
expected on quite general grounds, LDA locates the filled 
0-p states at too high and the empty Sr-d manifold at 
too low energies. The t2g manifold undergoes a strong 
quasi-particle renormalization with a concomitant shift 
of spectral weight to the lower Hubbard band, both of 
which are effects beyond the one-particle picture. 

Applying the GW approximation (see the spectral 
function in Fig. [T]d) increases the 0-p to Sr-d distance, 
placing both manifolds at energies nearly in agreement 
with experiment [5^. Most interestingly, however, a peak 
at 2.6 eV emerges from the d-manifold, which we find to 
be of Cg character. Indeed, the GW approximation en- 
hances the crystal field splitting and places the maximum 
of the Cg spectral weight at the location of the experimen- 
tally observed d^ d^ addition peak. 

Also displayed is the partial t2g contribution. Two 
interesting points are seen : the quasi-particle width is 
narrowed from the LDA value by about 0.5 eV and satel- 
lite structures appear below and above the Fermi level 
at ±3 eV, -|-5eV as well as at around 15eV. Within the 
(non-self-consistent) GW approximation, the self-energy 
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is related to ImW as follows 
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FIG. 1: (a) Orbital-resolved LDA density of states, 
in comparison to experimental spectra from photoemission 
(PES) [43, [131 and inverse photoemission (BIS)[43|. The lat- 
ter are reproduced also in panels (b-c). (b) Spectral function 
from GW (full orbital and t2g only) in comparison to exper- 
iments, (c) Spectral function from GW-I-DMFT (full orbital 
and t2g only) in comparison to experiments. The inset shows 
a low energy zoom. In the GW+DMFT spectrum, LHB de- 
notes the lower Hubbard band. See text for details. 
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Thus, a feature at luq in ImM^ causes a satellite feature 
at an energy ujq below or above the quasiparticle peak 
at Skn- In Fig. [2] we display the on-site t2g intra-orbital 
clement Wt2g of the fully screened interaction in the Wan- 
nier basis. Three notable features are clearly discernible: 
There is the usual high-energy electron-gas-like plasmon 
excitation at around 15 eV and, more remarkably, there 
is a strong excitation between 2-3 eV, as well as a shoul- 
der at '^4.5eV. Evidently, these energy positions directly 
correlate with the spectral features discussed above. It is 
important to note that the lowest feature in Wt2g mainly 
arises from a collective excitation within the t2g bands. 
This can be understood by comparison with the Hubbard 
interaction U{uj) of the low energy t2g system, as calcu- 
lated within the constrained RPA. In U{uj) the structure 
at 2~3 eV is absent, due to the elimination of the po- 
larization within the t2g bands. The shoulder at around 
4.5eV, however, partly persists, as its excitations also in- 
volve 0-2p states; we will come back to this point below. 

Solving the GW-f DMFT equations, even without the 
self-consistency cycle at the GW level, that is, for fixed 
and U{uj), is a hard task. We use the scheme 
recently introduced in Ref. 4]: the Green's function of 
the dynamical impurity model is obtained from a fac- 
torization ansatz G(r) = Gstat{T)B{T) where Gstat is 
the Green's function of a static impurity model with 
constant U=U{uj=0) 60], and the second factor B{t) = 
G{T)/Gstat{'T) is approximated by its value for vanishing 
bath hybridization A 4] : 



B{t) = exp 



with Kt(lo) 
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In the regime 
that we are interested in, namely when the plasma fre- 
quency is several times the bandwidth, this is an excellent 
approximation Q . 

Finally, we present our GW-(-DMFT results for the 
spectrum of SrVOa at inverse temperature j3=l{)eV~'^ : 
Fig.[Hc) displays the local projection, while Fig.|3]shows 
momentum dependent t2a spectra in comparison with 
ARPES measurements 39l |48| . The low-energy part 
of the spectrum is dominated by the t2g contribution, 
which, here, is profoundly modified with respect to pure 
GW results. A renormalized quasi-particle band dis- 
perses around the Fermi level : At the F point the peak 
is located at about -0.5 eV - a strong renormalization of 
the corresponding Kohn-Sham state which, at this mo- 
mentum, has an energy of -1 eV. At the X-point, the t2g 
bands are no longer degenerate, and surprisingly weakly 
renormalized xy/xz states are observed at 0.9 eV, while 
the yz band is located at nearly the same energy as at 



the r point, in agreement with ARPES. At binding en- 
ergies of 1.6 eV, ARPES witnesses a weakly dispersive 
Hubbard band, whose intensity varies significantly as a 
function of momentum [is*] . In the GW+DMFT spectral 
function the Hubbard band - absent in GW - is correctly 
observed at about -1.6 eV, and its k-dependent inten- 
sity variation (see Fig. [3]) is indeed quite strong. Previ- 
ous LDA-I-DMFT calculations placed the lower Hubbard 
band at larger negative energies (see e.g. [27]). This 
is owing to the fact that when using a static Hubbard 
interaction, a value of 4-6 eV, that is larger than the 
zero frequency limit of the ab intio U(uj=0)=3.6eV\l^. 
is needed to account for the observed transfers of spec- 
tral weight. LDA-|-[/(cl;)-|-DMFT calculations taking not 
only the ab initio value of the static component of U but 
also its full frequency dependence into account indeed 
reproduce the position of the lower Hubbard band quite 
well djliH. Analogously, in GW+DMFT^ the additional 
transfers from the dynamical screening jj, 0, HI] , yield a 
good description of the Hubbard band and the spectral 
weight reduction at the same time. 

At positive energies non-local self-energy effects are 
larger. Interestingly, our total spectral function, FigHKc), 
does not display a clearly separated Hubbard band. The 
reason is visible from the k-resolved spectra: the upper 
Hubbard band is located at around 2 eV, as expected 
from the location of the lower Hubbard band and the 
fact that their separation is roughly given by the zero- 
frequency value of U. The peak around 2.7 eV that ap- 
pears in the inverse photoemission spectrum [sol ] - com- 
monly interpreted as the upper Hubbard band of t2g char- 
acter in the DMFT literature - arises in fact from Cg 
states located in this energy range. The non-local self- 
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FIG. 3: Momentum-resolved t2g spectral functions from GW 
and GW-I-DMFT, compared to experimental angle resolved 
photoemission spectra of Refs. [3§] and [4^ . The GW spectra 
are calculated at zero temperature. 



energy effects lead, in the unoccupied part of the spec- 
trum, to overlapping features from different k-points and 
an overall smearing of the total spectral function. 

Akin to the correspondence between ImW and the GW 
spectrum, our approach to solve the GW-I-DMFT equa- 
tions allows for a transparent physical interpretation of 
the arising spectral properties. Indeed the spectral rep- 
resentation of the bosonic renormalization factor B{t) of 
Eq. [T] (displayed in the lower panel of Fig. [2|) is directly 
related to the density of screening modes ImU (u)) jixP' [4j . 
In this way, we can trace back the GW-I-DMFT satel- 
lite at -4.5 eV to the onset of "p-^ig excitations, discussed 
above for W and U . On the other hand, since the fea- 
ture below 3 eV in W is absent in U and -B, at least the 
local self-energy contribution to the corresponding peak 
within GW does not survive in GW-t-DMFT. The strong 
peak at 15 eV is the well-known plasma excitation, seen 
e.g. in electron energy loss spectra of SrTi03[53]. 

In conclusion, we have presented the first dynamical 
implementation of a combined GW and dynamical mean 
field scheme. Quite generally, materials with a "dou- 
ble LDA failure"- an inappropriate description of corre- 
lated states, and deficiencies of LDA for the more itin- 
erant states, such as an underestimated "pd-gap" - can 
be treated within our scheme. Application to a well- 
studied benchmark material, SrVOa, confirms the abil- 
ity of the approach to describe simultaneously Hubbard 
bands, higher energy satellite structures and corrected 
energy gaps, in an ab initio fashion. In the unoccupied 
part of the one-particle spectral function, we attribute 
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a strong peak evidenced in BIS spectra to the electron 
addition process into the eg states, in agreement with 
cluster studies 8]. We further predict the existence of 
a satellite feature at -4.5 eV, that may seem difficult to 
resolve experimentally because of its overlap with oxygen 
p-states. Nevertheless we would expect that a systematic 
PES study as a function of photon energy, due to the de- 
pendence of the cross section on orbital characters, could 
possibly resolve this feature. In fact, structures of similar 
origin may already have been observed in other materials, 
such as e.g. in VO2 54 1 and iron pnictides Q. 
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